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Abstract 

Starting with the many-body Schrodinger equation we derive new rate equa- 
tions for resonant transport in quantum dots hnked by ballistic channels with 
high density of states. The charging and the Pauli exclusion principle effects 



o 

CJ ' were taken into account. It is shown that the current in such a system dis- 



^f\ ' plays quantum coherence effects, even if the dots are away one from another. 

>■ 

A comparative analysis of quantum coherence effects in coupled and separated 
t^ ■ 

~j I dots is presented. The rate equations are extended for description of coherent 

O' 

t^^ ' and incoherent transport in arbitrary multi-dot systems. It is demonstrated 



that new rate equations constitute a generalization of the well-known optical 
Bloch equations. 



o 

O ■ I. INTRODUCTION 



Quantum transport in small tunneling structures (quantum dots) have attracted great 
attention due to the possibility of investigating single-electron effects in the electric current 
|1[]. Until now research has been mostly concentrated on single dots, but the rapid progress in 
microfabrication technology has made it possible the extension to coupled dot systems with 
aligned levels 0-0. An example of such a system is shown schematically in Fig. 1, where the 
coupled quantum dot (coupled wells) are connected with two separate reservoirs (the leads). 



In contrast with a single dot, the electron wave function inside a coupled dot structure 
is a superposition of electron states localized in each of the dots. As a result, the effects 
of quantum coherence would appear in electron current flowing through such a system. 
Usually these effects are treated in the framework of single electron approach. Although 
this approach is an appropriate one for coupled- well semiconductor heterostructures |^|-0, 
it cannot be applied for coupled dots. The reason is the dominant role of the Coulomb 
interaction that leads to the Hubbard-type Hamiltonian for a description of these systems 
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Fig. 1. Resonant transport through a coupled-dot structure. Tl^r denote the 
corresponding tunneling rate from (to) the left (right) reservoirs, and Oq is the 
interdot hopping amplitude. Ep denote the Fermi level in the left (right) 
reservoirs. Only those energy levels inside the dots that carry electron current 
are shown. 

The Coulomb blockade effects in electron transport through a single dot can be taken 
into account in the most simple way by using the "classical" rate equations [P|^T^. These 
equations describe the current through the dot in terms of balance between incoming and 
outgoing rates from (to) the leads. Basically the same classical approach has been used to 
calculate the conductance of coupled multi-dot systems [^,0. In these works the whole 
array of quantum dots has been treated as a single quantum system, where its many-body 
eigenstates were calculated by exact (numerical) diagonalization of the Hubbard Hamilto- 
nian. Then the leads (reservoirs) were incorporated through the rate equations for single 



dots [11]. Yet, such a procedure is correct only for small coupling with the leads. If the 



corresponding rates F are of the order of the interdot transitions fl (Fig. 1), they affect 
the diagonahzation of the Hubbard Hamihonian 0. Besides, the approach [^3|,|I^ is mainly 
numerical and the inelastic scattering within the array is accounted for phenomenologically. 
A different way to treat the quantum transport in multi-dot systems is to use the Bloch- 
type rate equations |T^|-P^ instead of the classical rate equations. In contrast with previous 
treatments ||l3| , p!^ this approach is valid also for strong coupling with the leads and it 



accounts the inelastic scattering effects inside the multi-dot system [|l^]. It proves to be 
much simpler than the other approaches, so that in many cases the analytical treatment 
of the problem is possible ||T6|JT7|1 . It is also important to note that this approach is not 
a phenomenological one: new generalized Bloch-type rate equations were derived from the 
many-body Schrodinger equation with the Hubbard Hamiltonian by integrating out the 



continuum states of the reservoirs 17 



An important advantage of the Bloch-type equations is a clear distinction between coher- 
ent ("quantum") and non-coherent ("classical") terms. The quantum coherence in the Bloch 
equations manifests itself in the non-diagonal density- matrix elements ("coherences"), cou- 
pled with the diagonal density- matrix elements ( "probabilities" ) |TP|. Such a coupling always 



takes place, whenever a carrier jumps between isolated states inside the system |]T^. Other- 
wise, the diagonal and non-diagonal matrix elements are decoupled and the Bloch equations 
turn into the classical rate equation. 

Actually Bloch equations describe quantum motion between isolated (non-orthogonal) 
states, which are directly coupled (as in the double- well potential). In general, these states 
can be separated by a medium with high density of states. An example of such a system is 
shown in Fig. 2, where two dots with aligned levels are separated by the ballistic channel. In 
this case a carrier transport from the left to the right reservoir proceeds through continuum 
states in the channel. Usually transitions from discrete to continuum states lead to dephas- 
ing, i.e. to destruction of quantum coherence. Yet, some coherent effects may survive. It 
is therefore important to establish of whether quantum coherence still affects the transport 
through separated dots, or more precisely: are the diagonal density-matrix elements still 



coupled with the non-diagonal matrix elements in the corresponding rate equations? 
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Fig. 2. Resonant transport through two quantum dots separated by a ballistic 
channel. Here Qi and fir denote the coupling of the left and right dots dot 
with the levels Ei and E^ in the left and in the right leads, ^m and Jlm denote 
the coupling of the left and the right dots with the level Em in the ballistic 
channel. 



This question leads us to a more general problem, of how to modify the Bloch rate 
equations for transport in coupled-dots for a general case of arbitrary distributed quantum 
dots. Similar to the Bloch equations one can anticipate that that the new rate equations 
would incorporate the classical and quantum descriptions in a very efficient way. Therefore 
these equations can provide most useful account of quantum coherence effects in different 
processes in a comparison with usual approaches [^ . 



In this paper we study the above mentioned problems taking as a basis the microscopic 
many-body Schrodinger equation with the Hubbard-type tunneling Hamiltonian, describing 
the entire system of the reservoirs and the quantum dots. The plan of the paper is the 
following. In Sect. 2 we describe the previously derived Bloch- type rate equations for 
coupled dots and some particular features of the coherent transport in these systems. In 
Sect. 3 we present the rate equations for separated dots and outline their derivation from the 
microscopic Schrodinger equation. The details of the derivation are described in Appendix. 
A comparison of coherent effects in quantum transport through coupled and separated dots 
is given in Sect. 4. The quantum rate equations for a general configuration of quantum dots 
are presented in Sect. 5. The last section is a summary. 



II. RATE EQUATIONS AND COHERENT EFFECTS IN COUPLED 

QUANTUM-DOTS 

We start with a review of the quantum rate equations for coupled multi-dot systems, con- 
nected with two reservoirs (the "emitter" and the "collector") and interacting with phonon 
reservoir. The entire system is described by the Hubbard-type tunneling Hamiltonian where 
the electron-electron interaction is taken into account by introducing the corresponding 
electrostatic charging energy. In case of large voltage bias one can reduce the many-body 
Schrodinger equation to the system of quantum rate equations by integrating out continuum 
reservoir states ^T^. As a result the following Bloch-type equations for the density matrix of 



the multi-dot system {cTaf^} are obtained, where a,f3 denote the isolated (non-orthogonal) 
states of the system: 

O'aa = "^ Z_/ ^a^l3i'^af3 " 0"/3a) " O'aa 2_/ rQ-»7 + 2_/ '^55^S~^a , (2.1a) 

/3 7 5 



&a/3 = i{Ef3 — Ea)(Jal3 + M /_/ (^a-f^j-^IS — z2 ^a-^S^rS/S 1 



~ 9*^0/3 z2 i^a^S + r/3^<5) + 2_/ '^-yS'^-yS-^a/S- (2.1b) 

Here fla^is is the amplitude of one-electron hopping that results in the transition between 
the states a and /3. The width Fq,^^ = 2Trp\Qa^^\'^ is the probability per unit time for the 
system to make a transition from the state \a) to the state I7) of the device due to the 
tunneling to (or from) the reservoirs, or due to interaction with the phonon bath, or any 
other interaction generated by a continuum state medium with the density of states p. Notice 
that Eq. ( p.la|) for diagonal elements has a form of classical rate equations, except for the 



first term. This term is generated by an electron hopping between the isolated levels, which 
results in the coupling with non-diagonal density-matrix elements. Therefore it is responsible 
for coherent quantum effects in electron transport. The nondiagonal matrix elements are 
described by Eq. ( p.lb|) . The last term in Eq. ( p.lbj ) appears only for systems with the 
number of isolated states participating in the transport is more than two. It describes the 



simultaneous conversion of the states 7 — ;> a and (5^/3, generated by the same one-electron 
decay to (from) continuum. 

The current flowing through the system is given by 

m=j:^..it)^R\ (2-2) 

where the sum is extended over states I7) in which the dot adjacent to the collector is 
occupied (we consider the electron charge e = 1). T^l' is the partial width of the state I7) 
due to tunneling to the collector (the right reservoir). 

It follows from Eqs. (|2.1| ) that the coherent effects do appear in the quantum transport 
whenever a carrier jumps from one to another isolated state inside the device. In the absence 
of such transition as, for instance, in resonant tunneling through a single dot, the diagonal 
and non-diagonal matrix elements are decoupled and the evolution of diagonal density-matrix 
elements is described by the classical rate equation. 

For an example we consider the quantum transport through a double-dot system at zero 
temperature [^, shown in Fig. 1. In order to diminish the number of equations we assume 
that the Coulomb repulsion does not allow for two electrons to occupy the same dot, i.e. 
Ei + Uii ^ Ep, where i = 1,2 and Ua is the corresponding charging energy. Yet, the 
interdot Coulomb repulsion Uu is much smaller, so the both levels Ei^2 can be occupied 
simultaneously. In this case there are four available states of the double-dot system (for 
simplicity we neglected the spin): \a) - the levels E12 are empty, \b) - the level Ei is 
occupied, \c) - the level E2 is occupied, \d) - the both level £'1^2 are occupied. Using 
Eqs. ( p.l|) with a,/3, ... = {a,b,c,d} we find the following equations describing the time 
evolution of the corresponding density-matrix elements 

C^aa = -^L(^aa + ^R(^cc , (2.3a) 

abb = ^L(^aa + ^'R(^dd + i^oicTbc - (Tcb) , (2.3b) 

CTcc = -^R(^cc - ^'l^cc - i^o(,(^bc " (^cb) , (2.3c) 

o'dd = -^'R(^dd + ^'l^cc , (2.3d) 
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&bc = i{E2 - Ei)abc + i^o{(ybb - cTcc) - ^{^'l + ^RJc^bc , (2.3e) 

which are supplemented with the probabihty conservation condition: J2i<^ii(t) = 1. Note 
that the interdot Coulomb repulsion results in a variation of the corresponding transition 
rate whenever both levels are occupied {T^ji is replaced by F^ ^, evaluated at the energy 
-^1,2 + Ui2 respectively |[T7|1). 



The dc current is / = Tjiaccit -^ oo) + T'^addit -^ oo), Eq. ( |2.2| ). In order to shorten the 
final expression we assume that the rates are weakly dependent on the energy, F'^ j^ = Tlr. 
Then solving Eqs. (|2.3|) in the limit t ^ oo we obtain 



Sl + rj e^TLTn/iTL + Tn)^ + ^l + Fz.FV4' ^ ' ^ 

where e = Ei — E2. 

It is interesting to compare Eqs. ( |2.3| ) with the corresponding classical rate equations 
where only transitions between the diagonal density-matrix elements take place. In this 
case the probabilities af,b and ace of finding an electron in the first and the second dot are 
coupled by the some rate Tm- One easily finds that Eqs. ( p.3|) become 

O-aa = -'^LO'aa + Fi?(Jcc , (2.5a) 

<^bb = FLCTaa + ^ R<^dd - Fj\/((T66 - CTcc) , (2.5b) 

O-cc = -Fi?acc - TlCTcc + ^MicThb - CTcc) , (2.5c) 

CTdd = -'^RC^dd + ^L<ycc ■ (2.5d) 

(Cf. with the analogues rate equations in [p2| , p3| ). Solving these classical rate equations in 
the limit t ^ oo we obtain for the (classical) dc current 



^l^bJ^m 



^L^R + ^L^M + r_RrM 



(2.6) 



Comparing Eq. ( |2.6| ) with its "quantum" counterpart, Eq. ( p.4| ), we find that the 
quantum-mechanical nature of the dc current displays itself in the Lorentzian-shape res- 
onance as a function of the levels disalignment, e = E2 — Ei. Such resonances were already 



observed in the coupled-dot systems 0. Yet, the existence of the resonance is not neces- 
sarily related to the quantum coherence effect. The later manifests itself in a peculiar and 
even a counter-intuitive dependence of the resonant current on coupling with the reservoirs. 
It follows from Eq. ( |2.4| ) that the current decreases when the coupling with the reservoirs 
increases. For instance, the peak value of the resonant current 1^=0 o^ ^l/^r ~^ for 
r^^ Tl- It implies that by increasing the penetrability of the barrier, connecting the sec- 
ond dot with the collector we diminish the total current. As an example we show in Fig. 3 
the peak value of the resonant current as a function of F/j for constant F^, and Qq (the solid 
line). This effect can be observed in coupled-dot structures by changing the corresponding 
gate voltage. 
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Fig. 3. Current through the coupled-dot structure, Fig. 1, as a function of 
Tr/Tl. The solid line corresponds to Eq. (2.4) for e = and r|^/4i7Q = 5. 
The dashed line is the result of the classical rate equations, Eq. (2.6), for 
Fm = Fl. 
For a comparison the corresponding "classical" current I^, Eq. (p^.6|), is shown in Fig. 
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by the dashed hne. 

The decrease of the resonant current shown in Fig. 3 is a result of destructive quantum 
interference during decay of a coherent superposition to continuum in the regime of strong 



couphng with the reservoirs p3| , p^ . A similar effect had been observed in electron decay 



from coupled quantum wells to continuum |^. It is quite interesting that such a counter- 



intuitive behavior of the resonant current looks as a manifestation of the Zeno effect p6 
It tells us that the quantum transitions between different states slow down when one of the 
states is continuously observed. Indeed, the increase of the right barrier penetrability leads 
to immediate tunneling of an electron to the collector whenever it arrives to the second dot. 
It can be considered as a continuous observation of the second dot state, which results in an 



effective electron localization in the first dot [27 



III. RATE EQUATIONS FOR TWO SEPARATED DOTS 

Now we are going to derive rate equations for quantum transport in mesoscopic systems 
with arbitrary configuration of quantum dots. As a generic example we consider quantum 
transport through two quantum dots separated by a ballistic channel. Fig. 2. The dots 
contain only isolated levels, whereas the density of states in the ballistic channel and in the 
emitter and the detector is very high (continuum). This system can be described by the 
tunneling Hamiltonian 

I m r i,j=l,2 

+ < Y^ ^ia\ai + Y^ ^mal^ai + ^ ^malcim + Y ^r-ajct2 + H.c. ^ (3.1) 

[^ I m m r ) 

Here the subscripts /, m and r enumerate correspondingly the levels in the left reservoir, in 
the (middle) ballistic channel and in the right reservoir. The spin degrees of freedom were 
omitted. In order to simplify the derivation we assumed that the intradot charging energy Ua 
is large, E12 + Ua ^ Ep. Thus only one electron can occupy each of the dots. However, the 
interdot charging energy U12 is much smaller, so it does not prevent simultaneous occupation 

9 



of the two dots. The same is assumed for the Coulomb repulsion between electrons inside 
the dots and the ballistic channel. Although we did not include this interaction in the 
Hamiltonian ( |3.1| ), it can be treated in the same way as the interdot interaction Ui2- As in 
the previous case we restrict ourselves to the zero temperature case, even though the results 
are valid for a finite temperature, as would be clear from the derivation. 

Let us assume that all the levels in the emitter, in the ballistic channel and in the 
collector are initially filled up to the Fermi energies Ep, Ep and Ep respectively. We call 
it as the "vacuum" state, |0). (In the following we consider the case of large bias, so that 
Ep ^ Ep,Ep). This vacuum state is unstable: the Hamiltonian Eq. ( p.l|) requires it to 
decay exponentially to continuum states having the form aja^lO) with an electron in the 
level El and a hole in the emitter continuum, aj^a^ |0) with an electron in the level Em in 
the ballistic channel and a hole in the emitter, and so on. The many-body wave function 
describing this system can be written in the occupation number representation as 



I^W) 



Lm 



M^) + J2 bii{t)a\ai + Y^ birait)al^ai + J2 hi{t)a20'i 

|0), (3.2) 



+ YK{t)alo'i + Yhi2w{t)a\a\aiav + Yl hiwr{t)a\alaiaii + . 

l,r KV l<l',r 

where b{t) are the time-dependent probability amplitudes to find the system in the cor- 
responding states described above. These amplitudes are obtained from the Shrodinger 
equation i|\E'(t)) = 7-^|\E'(t)), supplemented with the initial condition (&o(0) = 1, and all the 
other 6(0) 's being zeros). Using the amplitudes b{t) we can find the density-matrix of the 
quantum dots, a^j (t) , by tracing out the continuum states of the reservoirs and the ballis- 
tic channel. Here the subscript indices in a denote four states of the dots: i,j = {a, b, c, d}, 
where \a) - the levels Ei^2 are empty, \b) - the level Ei is occupied, \c) - the level E2 is 
occupied, \d) - the both level Ei^2 are occupied, and the superscript indices k,n denote the 
number of electrons accumulated in the ballistic channel and in the collector respectively at 
time t, Fig. 4. 
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Fig. 4. Electron states of the two separated dot structure, shown in Fig. 2. 
^L,R, '^M and Tm are the tunnehng rates between the dots and the reservoirs, 
and between the dots and the balhstic channel. The indices k and n denote the 
number of electrons penetrating to the ballistic channel and to the collector 
at time t. 
One finds 

l,m l<l',m,r 

<riTit) = j:\bn{t)\\ 4l''\t)= E \biu'Ut)\', 

I l<l',m 



(^bh (t) — E \blU'l"mr{t)\ , 



l<l'<l",m,r 



(1,0), 



<''\t) = j:biiimiit), aii''\t) = E buumu'mit) 



l<l',m 

(^bc (t) = E biu>in^,.{t)h*2iiq„mr{^) , 

l<l'<l",m,r 



(3.3a) 



(3.3b) 



(3.3c) 



The rate of electrons arriving to the collector determines the electron current in the 
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system. Therefore the current operator is / = i\H, iV^], where Nr = Y,r 'A-'^r i^ the operator 
for the total number of electrons accumulated in the right reservoir. Using Eqs. ( p.l|) , ( p. 21) , 
( p.3|) we find that the current I{t) flowing through the system is 

m = {^{t)\i\^>{t)) =Y.n [a^^-\t) + aS'")(t) + d(^")(t) + aff\t)\ (3.4) 

k,n 

As expected, I{t) is the time derivative of the total charge accumulated in the collector. 

It follows from Eq. ( p.4[ ) that the current /(t) flowing through this system is expressed 
in terms of the diagonal elements of the density-matrix cr(t). In order find the differential 
equations for a(t) we need to sum over the states of the reservoirs and the ballistic channel, 
Eqs. ( p.3|) . It can be done analytically by using the procedure similar to that in Refs. |1T7| , [18 



providing that the levels are not close to Fermi levels, \Ei — Ef\ ^ T, Fig. 4. The details of 
derivation are described in Appendix. As a result we obtain the following Bloch-type rate 
equations for the matrix elements of the density-sub matrix cr(t): 






^ cc 



. (fc 



+ 27rpMfiMl^M(o-L'~''"^ + a^r^'"^) (3.5a) 

- vrpM^M^Af (o-foc'"'' + o-^b'""*) (3.5b) 
-") = -(Fi + f M + F«)aL^'") + FMa£-^'") - 7rpMfiM^M(aL''"^ + a^'"^) (3.5c) 

""^ = -(Fm + Tm + rn)afr^ + T^a^^^ (3.5d) 

-^(^L + FM + fM + Ffi)aL'•"^ (3.5e) 

Here F denote the tunneling widths to the leads and the ballistic channel (F = 27rp|fip, 
where p is the corresponding density of states), and VLm, ^m are the hopping amplitudes 
between the left and the right dots and the states Em = -Ei,2 of the ballistic channel. 

Equations ( p.5|) have clear physical interpretation. Consider for instance Eq. ( p.5a|) for 



the probability rate of finding the system in the state a with k electrons in the ballistic 
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channel and n electrons in the right reservoir (Fig. 4a). This state decays with the rate T^ 
into the state h (Fig. 4b) whenever an electron enters the first dot from the left reservoir. 
This process is described by the first term in Eq. ( p.5a| ). On the other hand, the states 



h and c (Figs. 4b, 4c) with k — 1 electrons in the ballistic channel decay into the state a 
with k electrons in the ballistic channel. It takes place due to one-electron tunneling from 
the quantum dots into the ballistic channel with the rates Tm and Tm respectively. This 
process is described by the second and the third terms in Eq. (|3.5a| ). Also the state c (Fig. 
4c) with n — 1 electrons in the right reservoir can decay into the state a due to tunneling to 
the right reservoir with the rate Tr (the fourth term in Eq. (|3.5a|) ). The last term in this 



equation describes the decay of coherent superposition of the states h and c into the state 
a. It takes place due to single electron tunneling from the first and the second dots into the 
same state of the ballistic channel with the amplitudes VLm and VLm, respectively. Obviously, 
this process has no classical analogy, since classical particle cannot simultaneously occupy 
two dots. 

Equations ( p.5b| ), ( |3.5c[ ) and ( |3.5clD describe the probability rate of finding the system in 



the states where one of the dots or both dots are occupied. In the first case an electron can 
jump into unoccupied dot via continuum states of the ballistic channel. As a result, the states 
h and c can decay into linear superposition of the states h and c. This process is described 
by the last terms in Eqs. (|3.5b|) , ( |3.5c|) . Obviously, if the both dots are occupied, such a 



process cannot take a place. Therefore add is not coupled with the nondiagonal density- 
matrix elements, Eq. ( |3.5d|) . The last equation, ( p.5e|) describes the time-dependence of the 
nondiagonal density matrix element. It has the same interpretation as all previous equations. 
Equations (|3.5|) for the reduced density matrix cr(t) were derived starting from the wave 
function |\I'(t)), Eq. ( |3.2| ), instead of using the density matrix for the entire system. Of 
course, it makes no difference if the entire system is initially in the pure state, or all the levels 
of the reservoirs are occupied up to the corresponding Fermi energies. At finite temperature, 
however, the system is not initially in a pure state. Then one needs to perform the derivation 
using the Liouville (Landau-von Neumann) equation for the density-matrix. Yet, if energy 
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levels Ei^2 of the dots are far away from the Fermi-levels (T ^ Ep — -E'1,2, -E'1,2 — Ep), then 
the reservoir levels that carry the current {\Ei — -^1,2! ~ T) are deeply inside the Fermi sea. 
In this case one can consider these levels as fully occupied. Then neglecting the relaxation 
processes in the reservoir and assuming zero temperature inside the dots, we would arrive 
to the same Eqs. (|3.5|), derived from the pure state. 

Using Eqs. (p.5|) one finds for the total current, Eq. ( p. 41) 



/(t)=r^[a,,(t) + (Trfrf(t)], (3.6) 

where an = J2k,n ^u are the total "probabihties" . We can easily understand this result 
by taken into account that a^c + Udd is the total probability for occupation of the second 
dot and F^ is the rate of electron transitions from this dot to the (adjacent) right reservoir, 
Eq. (13). 

In order to find differential equations for aij we sum over k,n in Eqs. (|3.5| ). Then 
we obtain the following Bloch-type equations, which describe the time-dependence of the 
density-matrix for separated dots 

(^aa = -^L(^aa + ^ MCTbb + ^MCTcc + ^R(^cc + '2T<'PM^M^M{(^bc + (^cb) (3.7a) 

0"fe6 = -^M<^bb + ^L<yaa + ^M<^dd + ^R<^dd - T^PM^M^Mi'^bc + (^cb) (3.7b) 

C^cc = -(Xl + f A/ + Tji)acc + ^MfTdd - T^PM^M^M{.(^bc + (^cb) (3.7c) 

Odd = -C^M + f A/ + ^R)(ydd + ^L<ycc (3.7d) 

dfoc = i{E2 - Ei)abc - TtpM^M^Micfbb + (^cc) + 27rpMQM^M<ydd - -^^totCTbc, (3.7e) 

where Ttot = El + Tm + Tm + ^r- 

The stationary (dc) current / = I{t -^ oo) Eq. ( |3.6| ) can be easily obtained from Eqs. ( p.7|) 
by taken into account that a^ -^ for t ^ oo. Then Eqs. ( |3.7| ) turn into a system of linear 
algebraic equations, supplemented by a probability conservation condition (Xaa + cTfeb + <^cc + 
add = 1- Consider, for example, the case of the same partial widths of the levels -Ei,2, i-e. 
Tl = Tji = Fo and Tm = f m- The latter implies 2ti p^VLm^m = ^Tm, since the amplitudes 
Qm, ^m can be of the opposite signs. Solving Eqs. (|3^ ), ( p.7|) one finds for dc current 
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2e2(ro + Tm) + 2To{Tl + 3ToTm + STl,) ' ^ " ^ 

where e = Ei — E2. 

Similar to the coupled-dot case, Eq. ( p. 41 ), dc current in separated dots displays the 
Lorentzian shape resonance as a function of e and the same peculiar dependence on the 
coupling with the leads. The latter manifests the quantum-coherence effects in separated 
dot systems. We discuss these effects in the next section. 

IV. COHERENT EFFECTS IN TWO-DOT SYSTEMS 



It follows from the rate equations (|3.7|) that the diagonal density-matrix elements are 
coupled with the non-diagonal terms {ate, CTcb), similar to Eqs. ( |2.3| ) for a coupled-dot system. 
It means that the coherent effects should survive even in separated quantum dots. Indeed, 
let us consider the dependence of the resonant current on a coupling with the reservoirs. As 
an example, we show in Fig. 5a the peak value of the resonant current (e = 0) obtained from 
Eqs. (p.6|), ( |3.7| ) as a function of Tr at fixed Tl = Tm- One finds that the current through 
the dots separated by a ballistic channel displays the same peculiar behavior with Tji, as in 
the case of coupled dots. Fig. 3. Namely, the current decreases when the penetrability of 
the right barrier increases. Such a behavior of the resonant current is a clear manifestation 



of the quantum coherence effects, as explained in Sect. II (see also ||23|-p5|). It also can be 
verified experimentally by changing the voltage on the gate, connecting the right dot with 
the collector. 

On the other hand, the dependence of the resonant current on the width Tm is quite 
different. Fig. 5b. One finds that the current increases with Tm- It is not surprising, since 
Tm plays a role of coupling between the two dots. In fact, a similar behavior with the 
coupling ^0 displays the resonant current in the case of coupled dots, Eq. (|2.4[). 

Despite the analogues quantum coherence effects, the rate equations ( |3.7| ) are different 
from their counterpart for the coupled-dots, Eqs. (|2.3|). One finds that the diagonal matrix 



elements in Eqs. (|3.7|) are coupled with the real part of the non-diagonal matrix element, 

15 



i8=qA; 



^- 1 




(a) 



^e= 


U^^L 








0.15 






0.125 
0.1 






rR=ri 




0.075 










. 05 










0. 025 












5 


10 


15 


20 



(b) 



rMAl 



Fig. 5. DC current through the separated dot structure, Fig. 2, as a function 
of F/j/r^ for Fjv/ = F^ (a), and as a function of Fjv//F/^ for F/j = F^ (b). 

whereas the corresponding terms in Eqs. (|2.3| ), are coupled with the imaginary part of the 
non-diagonal matrix element. Notice also that the coupling is determined by the hopping 
amplitude fi in Eqs. ( p.3| ), while the corresponding coupling in Eqs. ( p.7| ) is proportional to 

It is important to point out that the coupling with the non-diagonal density-matrix 
elements in Eqs. (|3.71) does not decrease with the separation distance between the dots, 
although the hopping amplitudes Qm, ^m do decrease. Indeed, using semi-classical expres- 
sions one finds \Qni\ = (I/^tiTm) exp(— S') and \Q!^i\ = {l/^/^Y^M) exp(— S), where S{S) is 
the action under the barrier, separating the first (second) dot from the ballistic channel. It 
can be written as 5 = J^.^ \p{x)\dx, where \p{x)\ - 



'2m[V{x) — E] and Xjj, are the classical 
turning points, ri 2 are the classical periods of motion in the first (second) quantum well 
and tm is the classical periods of motion in the ballistic channel. It implies that tm — ^ cx3 
when the length of the ballistic channel increases. Yet, tm is canceled out from the cou- 
pling between the diagonal and non-diagonal matrix elements, since the ballistic channel is 
effectively one-dimensional so that the density of states is pm = TA//27r. Finally we obtain 

1 



27ipM^M^ 



M 



'2-K^TYfi 



exp[-(5 + 5)]. 



(4.1) 



16 



Therefore, the quantum coherence effects in dc current could survive even for long ballistic 
channels. Obviously, it is true only for ideal one- dimensional channels, where all the flux 
coming from the first dot arrives the second dot. It should be also pointed out that electron- 
electron scattering inside the ballistic channel (not taken into account in these calculations) 
may reduce the interference effects. Yet, we expect that such a process would not diminish 
the coupling (|4.1|) , but increase the "dissipation" width Ttot in Eq. (|3.7e| ). This problem, 
however, needs a special investigation. 

It is interesting to make a comparison between peak values of the dc current in coupled 
and separated dots. In both cases the current reaches its maximal value for Tl = T^ = Fq 
and e = 0, Eqs. ( p.4|) , (pl8|). One finds for the coupled dots 



Fo fin ^0 



^, for fio>ro, (4.2) 



""''' 2 QI + F2/4 2 
while for the separated dots 

Imax = -7ror2 , or T — ZT2 ~^ "«"' ^°^ ^A-f ^ ^0, (4.3) 

^ -jJ- Af + -JJ- O-L M + J- 

Thus the peak value of dc current in the coupled dots is three times larger than that in the 
separated dots. 

It is quite natural to associate quantum coherence effects in double-well systems with 
quantum oscillations of an electron between the wells. Indeed, an electron in a coupled-dot 
system, detached from the emitter and the collector, oscillates between the states (a) and 
(b). Fig. 6. These oscillations are reproduced by the same rate equations that describe 
quantum transport in coupled-dots, Eqs. ( p.3|) , but now with T^ = F/j = 0. We then obtain 
for the probability of finding an electron in the states (a) and (b), daait) and <7bb(t), Fig.6: 

aaa = i^o{<yab " (^ba) (4.4a) 

&bb = -i^oiaab - CTba) (4.4b) 

aab = i{E2 - Ei)aab + ^fio(o"aa " (^bb)- (4.4c) 
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Fig. 6. Electron states of the coupled-dot structure of Fig. 1, isolated from 
the leads. 
Solving this equations with the initial condition o"aa(0) = 1 and cr(,b(0) = 0, we obtain 



(^aa{t) 



nl + eV4 



(4.5) 



where cui = (l/2)w4r2g + e^. As expected, an electron oscillates between the dots with an 
amplitude 1/(1 + e^/^^l). 

The situation is different, however, if the dots are separated by a ballistic channel. Fig. 7. 
Indeed, the corresponding rate equations for the probabilities aaaif) and o"b(,(t) of finding an 
electron in the states (a) and (b) of Fig. 7 are obtained from Eqs. ( P^ ) for F^ = F^ = 0. 
Consider for simplicity the case of ^2^ = ^m- Then Eqs. ( ^.71 ) become 



1 
O'bb = — r^/CTfefe — -FM(0"afe + CTfja) 

<^ab = i{E2 - Ei)aab - -Fa/ (0"aa + (^bb) - ^M<yab 



(4.6a) 
(4.6b) 
(4.6c) 



Solving these equations with the initial condition craa(O) = 1 and ^^^(O) = 0, we obtain 

(Taait) = '' ^ '; exp(-FMt), (4.7) 

J- M ^ 



where UJ2 = (1/2)yF|j — e^. It follows from Eq (f4.7|) that <7aa(t) does not display any 
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Fig. 7. Electron states of the two separate dot structure of Fig. 2, isolated 

from the leads. 
oscillations. As an example, we show in Fig. 8 the probability of finding an electron in 
the left dot as a function of time in the case of coupled and separated dots (Figs. 6, 7) for 
e = Tm = f^o- The solid line corresponds to Eq. ( [4.51 ), and the dashed line to Eq. ( [4.7] ). No- 
tice that aaa = o'bb -^ 1/4 for t ^ oo when the levels are aligned (e = 0), and aaa = crtb ^ 
for t ^ oo when e ^ 0. It means that an electron decays into the ballistic channel. 



Oaa(t) 




Fig. 8. The occupation of the left dot in Figs. 6 and 7 as a function of time 
for e = Tm = ^o- The solid line corresponds to coupled dots, and the dashed 
line to separated dots. 
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The reason for a dissipative behavior of cjaa in the case separated dots is a couphng of 
the dots with continuum states of the balhstic channel. The latter leads to dephasing that 
destroys quantum oscillations. It is thus rather remarkable that the quantum coherence 
effects would "reappear" when the same system is connected with the leads, Fig. 2. 

V. GENERAL CASE 

The rate equations (|3.7|) , describing electron transport in separated dots can be extended 
to any multi-dot system. By applying the same technique of integrating out the reservoir 



states as in Ref. ^j and in Appendix, we arrive to the rate equations for the density-matrix 
(Jaj3 of the multi-dot system. These equations can be written as 

V 7 7 / 

~ Z^ '^ Pi'^ay^j^S^S^IS + CTj/S^j^S^S-^a) +Z^'^P {^j^a^S^fS + ^j^l3^5^a)<^jS ■ (5.1) 
7,5 7,<5 

Here |a), |/3), . . . are the states of the multi-dot system in the occupation number represen- 
tation, and ^a^f3 denotes one-electron hopping amplitude that generates a — > /5-transition. 
We distinguish between the amplitudes Cl and Q of one-electron hopping among isolated 
states and among isolated and continuum states, respectively. The latter transitions are of 
the second order in the hopping amplitude ~ Q^. These transition are produced by two 
consecutive hoppings of an electron with the same energy across continuum states with the 
density of states p. The symbol J2' means that only those transitions '-/—>■ 6 ^ {a, (3) are ac- 
counted for in the sum, where the number of electrons in continuum states (in the reservoirs 
and ballistic channels) stays the same in the initial and final states (c.f. Eqs. ( |3.5| )). 

It is rather easy to verify that Eqs. { f).l\ j coincide with Eqs. ( p.7| ) for «,/?,... = {a, b, c, d}, 
which are the states of the separated dot system, shown in Fig. 5. Let us compare now 
Eqs. ( |5.1|) with the Bloch rate equations (|2.1| ) for quantum transport in coupled dots. We 
begin with the equations for diagonal density-matrix elements, (3 = a. The first term in 
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Eq. ( |5.1| ) is zero, since Ea = Ep. The first term in Eq. (|2.1a|) and the second term in 
Eq. ( [5.1| ) have the same form and describe the coupling with non-diagonal matrix elements 
generated by one-electron hopping between isolated states. The remaining terms in both 
equations look differently. However, one can easily realize that in the case of coupled dots 
the only possible transitions for diagonal matrix elements are those corresponding to 7 = a. 
Then the third term of Eq. (|5.1| ) becomes 



^ 2'KpVLa^sVLs^a(yaa = (^aa ^ '^a^6- (5.2) 

<5 <5 



It coincides with the second ( "dissipative" ) term of the Bloch rate equation ( p.la|) . The last 



term in Eq. ( |5.1| ) describes the reverse process, i.e. all possible conversions of the states 7, 
5 into the state a. Again, this term turns into the third term of Eq. ( p.la| ) for 7 = 5. 



The same relation can be traced for the non-diagonal density-matrix elements. We find 
that the third term of Eq. ( |5.1| ), which describe all possible decays of the states a and (3, 
reproduce the third term in Eq. (p.lb|) for 7 = /? and 5 = a, respectively. The fourth term 



in Eq. (|5.1|) describes the reverse process, i.e. the conversion of a^s -^ cTa/s produced by two 
consecutive hoppings through continuum state media. In particular, if Q-y^a = ^s^p, these 
terms reproduce the last term of Eq. (|2.1b|) , where T^s->ai3 = 2'KpQ.y^a^s->i3- 



VI. SUMMARY 

In this paper we derived quantum rate equations, which provide the most simple and 
transparent way for a description of the both coherent and incoherent electron transport in 
quantum dot systems. In the beginning we considered a system of two quantum dots linked 
by a ballistic channel and connected with the emitter and collector reservoirs. Starting 
with the many-particle wave function and integrating out the continuum states, we have 
obtained the equations of motion for the density submatrix of the two-dot system in the 
occupation number representation. In spite of the quantum dots might be away one from 
another, we found that the diagonal density-matrix elements are still coupled with the non- 
diagonal density-matrix elements, similar to Bloch equations for double-well systems. The 
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essential difference, however, is that the diagonal matrix elements in the Bloch equations are 
coupled with the imaginary part of the non-diagonal density-matrix elements, while for the 
separated wells the corresponding terms are coupled with the real part of the non-diagonal 
density-matrix elements. 

As expected, the coupling between diagonal and non-diagonal matrix elements gener- 
ated quantum coherent effects in electron transport. For instance, we found that due to 
destructive quantum interference the dc current in separated dots is reduced when the bar- 
rier penetrability increases. This effect would survive even if the dots are largely separated, 
providing that they are linked by ideal one-dimensional ballistic channels. 

In spite of the quantum coherence effects found in dc current, the same system of sep- 
arated dot, detached from the emitter and collector reservoirs does not show any quantum 
oscillations. It is quite different from a coupled-dot system with aligned levels, in which an 
electron oscillates between the dots. 

The rate equations for two separated quantum dots are extended to a general case of 
multi-dot system, in which the dots are either directly coupled, or interconnected via ballistic 
channels. These new rate equations generalize the well-known Bloch equations, describing 
time-evolution of the density-matrix of coupled multi-well systems in the presence of a 
dissipative media. We thus expect that the applicability of our generalized rate equations is 
not restricted by quantum transport in multi-dot system only, but these equations would be 
very useful for various physical problems where quantum coherence and classical dissipation 
effects do interplay. 
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APPENDIX A: DERIVATION OF RATE-EQUATIONS FOR SEPARATED DOTS 



We demonstrate of how the sum over continuum states in the density matrix, Eqs. 
can be performed analytically leading to the rate equations for two-dot systems. In order 
to simplify the algebraic transformations we use the Laplace transform for the Shrodinger 
equation i|\E'(t)) = 7Y|\E'(t)). Then the amplitudes b(t) in the wave function (|3^ ) are replaced 
by their Laplace transform 

b{E) = / e'^%{t)dt , (Al) 

Jo 

Substituting Eq. ( pJ.2D into the Shrodinger equation we obtain an infinite set of coupled 

equations for the amplitudes h{E): 

Ebo{E)-J2nME)=t (A2a) 

(E + Ei- E{)bn{E) - nME) - E ^mbim{E) = (A2b) 

m 

{E + Ei- EjbUE) - nJu{E) - nj2i{E) -Y^VlA^u'iE) = (A2c) 

I' 

{E + Ei- E2)b2i{E) - J2 ^Jim{E) - E ^rbir{E) - ^ fi,6m'(^) = (A2d) 

m r l> 

{E + Ei + El, -Ei~E2- Ui2)bi2ii'{E) - nib2i{E) + VtihuiE) - 

- E ^mblll'miE) - J2 ^mhu'miE) - ^ ^rbm'r{E) = (A2e) 



Notice that due to the Pauli principle an electron can return back only into unoccupied 
states of the emitter. As a result, the summation over the emitter states does not appear in 
the corresponding terms of Eqs. (|A2|) (in the second term of Eq. ( |A2b|) , in the second and 
the third terms of Eq. (|A2c|) , and so on). 



Eqs. (|A2|) can be substantially simplified. Let us replace the amplitude b in the term 



J2 ^b of each of the equations by its expression obtained from the subsequent equation. For 
example, substitute bii{E) from Eq. ([A2b|) into Eq. ( |A2a]) . We obtain 



^_y- ^liEi 



Y E + Ei-Ei 



bo{E)-J:^f§^bUE) = ^, (A3) 
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where Qi = Ql{Ei) and Qm = ^M{Em). Since the states in the reservoirs are very dense 
(continuum), one can replace the sums over / and m by integrals, for instance J2i ^^ 
J Pl{Ei) dEi , where Pl{Ei) is the density of states in the emitter. Then the first sum in 
Eq. ( [A^ ) becomes an integral which can be split into a sum of the singular and principal value 
parts. The singular part yields -iQ{EJ^ + E - Ei) Tl/2, where Tl = 27rpi(Ei)|fii(Ei)|2 is 
the level Ei partial width due to coupling to the emitter. Let us assume that Ep ^ Ei ^ 
Ep , i.e. the energy level is deeply inside the band. In this case the integration over Ek^^-^- 
variables can be extended to ±00. As a result, the theta-function can be replaced by one, 
and the principal part is merely included into redefinition of the energy Ei. The second sum 
(integral) in Eq. (^) proves to be negligible small. Indeed, let us replace him -^ b{Ei, E^, E), 
and assume weak energy dependence of Q on Ei. Then one finds from Eqs. ( [X^ ) that the 
poles of the integrand in the E'j-variable are on one side of the integration contour, and 
therefore this term vanishes. In general, any terms of the type / ■ ■ ■ dEs ■ ■ ■h{- ■ ■ ,Es,- ■ ■){E + 
■ ■ -^Es)^^ -^ 0, whenever the integration over the i^^^-variable can be extended to ±00. We 
shall imply this property in all subsequent derivations. Notice that these results are valid 
also for non-zero temperature, providing that T ^ Ep — -E'1.2, -^1,2 ~ Ep^ . 

Now we apply analogous considerations to the other equations of the system (^). It 
follows from Eqs. ( [A^ ) that the Coulomb repulsion just shifts the energy levels in Eq. ( |A3D , 
-^1,2 -^ -E'1,2 + Uij. If the shifted level is outside the band, the singular part of the integral 
vanishes and therefore the corresponding width F is replaced by zero. If, however, the 
shifted energy stands deeply inside the band, all the previous treatment remains the same. 
A problem appears only if the shifted energy in near the Fermi level. This case is not 
considered here. Finally we arrive at the following set of equations: 

{E + zFi/2)6o = i (A4a) 

{E + Ei-E^ + iTm/2)K - Qik + tTcpM^tM^lMhi = (A4b) 

{E + Ei-Em + iTL/2)bim - ^mbii - ^J>2i = (A4c) 

{E + E1-E2 + iVj2 + lYull + i^RJiyoii + mpM^M^ubu = (A4d) 
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{E + Ei + Ev -Ei-E2- Uu + ^r;,,/2 + tT\,/2 + tTy2)~h2U' - ^uhi + ^ihv = 0, (A4e) 



where Vm = '2'TiPai{Ei)\Qm{Ei)\'^ and Tm = 2'kpm{E2)\CIm{E2)\'^ are the partial widths of 
the levels Ei and E2 respectively due to coupling to the ballistic channel and pm is the 
density of states in the ballistic channel. Tn = 27rpij(£'2)|^i?(-E'2)P is the level E2 partial 
width due to coupling to the collector. The Coulomb interdot repulsion Uu just modifies the 
corresponding width (F is replaced by F', evaluated at the energy £'12 + f/12) whenever both 
dots are occupied |T^. A similar modification would take place if we include the Coulomb 
repulsion between electrons inside the dots and the ballistic channel. In the following we 
assume that widths are weakly dependent on the energy, so F' = F everywhere. 

The density matrix elements, Eqs. ( p.3|) , are directly related to the amplitudes b{E) 
through the inverse Laplace transform 

l...,m,...,r... 



Using this equation one can transform Eqs. ([A4|) for the amplitudes b{E) into differen- 



tial equations for the probabilities a^^''^\t). Consider, for instance, the term a^^' (t) = 
Y.I \bii{t)\'^, Eq. ( g^Sbl) . Multiplying Eq. ( |A4bD by h\i{E') and then subtracting the complex 
conjugated equation with the interchange E ^^ E' we obtain 
dEdE' 



47r2 



Y^{{E'-E- trM)bu{E)bUE') - Qi[hi{EME') - bUE%{E)] 



- iTipM^M^M[bii{E)bUE') + bUE')b2i{E)]}e'^'''~'^'^' = 0. (A6) 
Substituting 

''^ ^~ E + Ei-Ei + tTM/2 ^ ^ 

into Eq. ( |A6|) we can carry out the E, £"-integrations thus obtaining 

4r = -TMaT^ + F.air^ - np^nM^MiaS^''' + aT\ (A8) 



It corresponds to Eq. ( ^.Saj ) for k = n = 0. Applying the same procedure to each of equations 



( |A4c| ) we obtain the Bloch-type equations, (|3.5| ) for the density matrix element o"*^^'"'^(t). 
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